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Abstract. In this paper, by explicitly calculating the principal symbols of pseudodif- 
ferential operators and by applying Homander's spectral function theorem, we obtain 
the Weyl-type asymptotic formulas with sharp remainder estimates for the counting 
functions of the two classes of biharmonic Steklov eigenvalues Aj. and fif^ in a smooth 
bounded domain of a Riemannian manifold. This solves a longstanding challenging 
problem. 



1. Introduction 

Let {A4 , g) be a C°° Riemannian manifold of dimension n with a positive definite metric 
tensor g, and let 17 C be a bounded domain with C°° boundary dfl. Assume g is a. 
non- negative bounded function defined on dfl. We consider the following two biharmonic 
Steklov eigenvalue problems: 

( A^u = in n, 

(1.1) } u = on on, 

[ Agu + Apff = on 

and 

r A^v = in n, 

(1.2) 1^ = on an, 

[ 9(A^_^3^3^^o on 9^!, 

where v denotes the inward unit normal vector to dCl, and Ag is the Laplace-Beltrami 
operator defined in local coordinates by the expression. 

Here \g\ := det{gij) is the determinant of the metric tensor, and g*^ are the components of 
the inverse of the metric tensor g. 
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(1.1) and (1.2) are the biharmonic Steklov eigenvalue problems (see [6], [10], [17], [19], 
[26] . [4Tj and 03]). In each of the two cases, the spectrum is discrete and we arrange the 
eigenvalues in non-decreasing order (repeated according to multiplicity) 

< Ai < A2 < • • • < Afc < • • • , 

= Ml < M2 < • • • < Mfc < • • • • 

The corresponding eigenfunctions on dfl are expressed as ^7 , ^7 , • • • , , ■ ■ ■ ; vi, V2, ■ ■ ■ , 
Vk, - ■ ■ ■ It is clear that Afc and fik can be characterized variationally as 

Afc ^ ' ^ ' , = max inf ' ' , fc = 2,3,4,--- 

JaVl^\ du ) codim(j^)=k-i 0#§^ei2(sn) Jdn'^\di^l 



and 



\ = max inf , fc = 2,3,4,--- 



where H"^{rt) is the Sobolev space, and where dx and ds are the Riemannian elements of 
volume and area on and 951, respectively. 

The boundary value problems (|l.ip and (jl.2p have very interesting interpretations in 
theory of elasticity. We refer the reader to [TU], [IS] and [13] for more details. In view 
of the important applications, one is interested in finding the asymptotic formulas for Afc 
and /ife as fc — > 00. Let us introduce the counting functions A{t) and B{t) defined as the 
numbers of eigenvalues Afc and /ifc less than or equal to a given r, respectively. Then our 
asymptotic problems for the eigenvalues are reformulated as the study of the asymptotic 
behavior of A{t) and B{t) as t — > +00. 

The simpler harmonic Steklov problem was first introduced by V. A. Steklov for bounded 
domains in the plane in [36 . This problem is to find function v satisfying 



(1.3) 



AgW = in ri, 

§ + i^gv = on on, 

where 77 is a real number (The function v represents the steady state temperature on such 
that the flux on the boundary is proportional to the temperature). The harmonic Steklov 
spectrum of the domain is also called as the spectrum of the Dirichlet-to-Neumann map 
(see [7], [11]) or [37]. For the harmonic Steklov eigenvalue problem ()1.3p . in 1955 Sandgren 
[28] established the asymptotic formula of the counting function N{t) = if{k\rik < t}: 

(1.4) N{t) = , / g-^-^ds + o(t"^1) as r ^ +00, 



{2n) 



where a;„_i is the volume of the unit ball of M" ^. In the case that Riemannian manifold 
Ai and the boundary of f2 are smooth, the author |21j further gave a sharp remainder 
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estimate for the counting function of the harmonic Steklov eigenvalues {'r]k}'^^i- 

N{r) = j^^ g''-\x)dx + 0{t"-') as t ^ 

For the biharmonic Steklov eigenvalue problem (jl.ip with general domain, in [2T\ the author 
also established the leading asymptotic formula with remainder o(r"~^) as t — > +00. 

For the Dirichlet and Neumann eigenvalues of the Laplacian on a bounded domain f2 C 
R", H. Weyl ([44], [45]) in 1912 proved the following asymptotic formula which answered 
a question posed in 1908 by the physicist Lorentz: 

(1.5) 7Vzp(r) = (27r)~"w„(vol(17))T"/2 j_ o(t-"/2) as t +00, 

where N^{t) = if{k e N|afc < r}, N+{t) = #{fc e N|/3fc < r}, and < ai < 02 < • • • < 
ak < • • • and — f3i < (32 < ■ ■ ■ < (3k < ■ ■ ■ are all the Dirichlet and Neumann eigenvalues 
on n, respectively. As far back as in 1912, H. Weyl [17] conjectured (see also, Clark |S]) 
that the second term of the asymptotic formula for N^^ir) contain an (n — l)-dimensional 
measure ('area') of the boundary dfl, i.e., 

(1.6) N^ir) = (27r)-"w„(vol(17))T"/2 T J(27r)-"+ic^„_i(vol(af^))T("-i)/2 

+o(t("-i)/2) as t +00. 

In 1980, Ivrii 1201 proved this conjecture for domains having smooth boundary under the 
following condition regarding the billiard trajectories of fl, where the billiard trajectory in 
is taken with the usual reflections at the boundary (see also p. 100 of 0). Melrose [25] 
independently obtained the second term in Weyl's conjecture for manifolds with concave 
boundary. 

Note that it was already observed by Avakumovic [5] that for the Laplacian on the 
sphere §" , the high multiplicities of the eigenvalues make it impossible to improve (|1.5p to 
(|1.6p . Seeley (see, [30] and [31]) in 1980 gave the sharp asymptotic formula of the counting 
function (see also [1], [5], [37]): 

(1.7) iV=p(r) = (27r)-"a;„(vol(f7))r"/2 + 0{t^"-^'>/^) as t ^ +00. 

Applying the sharp asymptotic result ()1.7p . Sogge invented the unit band spectral projec- 
tion operator (see [33] and [34]) and established the well-known asymptotic estimates of 
the mapping norm Hxt-Hl^-^Lp (2 < p < 00) (cf. [31]). Howerver, it has been a longstanding 
challenging problem to get the sharp Weyl-type asymptotic formulas for the biharmonic 
Steklov eigenvalues (see [21]). 

In this paper, by explicitly calculating the principal symbols of the corresponding pseu- 
dodifferential operators for the problems (|l.ip and (|1.2p in the boundary of a C°° bounded 
domain, we obtain the sharp asymptotic formulas for the counting functions A(t) and 
B{t), respectively. The main results are the following: 

Theorem 1.1. Let {Ai, g) he an n- dimensional C°° Riemannian manifold, and letQ d A4 
be a bounded domain with C°° boundary dfl. Then 

, . ^n— 1 r 

(1.8) A{r) = J^^ g'-'ds + 0{r"~') as r ^ 
Moreover, the above remainder estimate is sharp. 
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Theorem 1.2. Let {M,g) be an n- dimensional C°° Riemannian manifold, and letQ, (Z Ai 
be a bounded domain with smooth boundary dfl. Then 

(1.9) Bir) = / ^""'^^ + ^(^""') «^ ^ ^ +°°- 

(VIdtt)" ^ JoD 

Moreover, the above remainder estimate is also sharp. 

The plan of the paper is as follows. In Section 2 we give some definitions and lemmas. 
In Section 3, by a key technique we calculate the principal symbols of the corresponding 
"Neumann-to-Laplacian map" and "Dirichlet-to-Laplacian derivative map" . Section 4 is 
devoted to the proofs of the sharp Weyl-type asymptotic formulas for A{t) and B{t). In 
Section 5, we give two counterexamples, which show that Theorem 1.1 and 1.2 cannot be 
improved. 



2. Some definitions and lemmas 



Definition 2.1. If U is an open subset o/M", we denote by S"" = S""([/,R") the set of 
all p G C°°(i7, M") such that for every compact set K <Z U we have 

(2.1) \DiiD1p{x,i)\<CK.,o.A^ + nr-\''\, xeif, eeK" 

for all a, (3 (z N" . The elements of S™ are called symbols of order m. 

It is clear that S*™ is a Frechet space with semi-norms given by the smallest constants 
which can be used in (j2.1|) (i.e., 



(2.2) 



\\P\\K,a.,P 



sup 



(i?fi?^"p(x,o)(i + iei)i"i"" 



Let p{x,£,) € S*™. A pseudo-differential operator in an open set U C 
defined by a Fourier integral operator (cf. [13,): 

1 



is essentially 



(2.3) 



P{x,D)u{x) 



(27r) 



p(x,e)e^<^'«>^(e)de 



Here u e C^{U) and u{^) = Jjg„ e ^^^'^'^ u{y)dy is the Fourier transform of u. 

Definition 2.2. A pseudodifferential operator P with its symbol p in S*™ is called classical 
or poly homogeneous if there is a sequence of symbols pj G 5'™^-', j = 0, 1, 2, • • • , such that 
Pj{x,tCj = t"^-^pj{x,^) for t > 1, \l\ > 1, and 



N 



d^,di[p{x,0-Y.pA^^^) 



< Ca,l3,N\i\' 



i-N-l-\a\ 



for all a,/3, \^\ > 1 and all integers N > 0. In this case the notation p ^ ^^oPj used. 
The function po is known as the principal symbol of pseudodifferential operator P, and the 
class of such symbol is denoted by S™ . 

Given a diffeomorphism t : t/i — > J72, from one open set Ui C M" onto another open 
set U2 C M", the induced transformation t* : C^{U2) C^{Ui), taking a function u to 
the function not, is an isomorphism and transforms C^{U2) into C^{Ui). Let Pi he a 
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pseudodifFerential operator on Ui and define P2 : C§°{U2) Coo (£^2) with the help of the 
commutative diagram 

-I I'- 

C^{U2) C^{U2) 

i.e., 

(2.4) P2U ^ [Pi{u o l)] o . 
(|2.4p can also be written as 

P2U= {l-^)*Pi{i*u). 

It follows from this that P2 is also a pseudodifferential operator on 1/2- 

Let be a smooth n-dimensional Riemannian manifold (of class C°°). We will denote 
by C°°{M) and C^{M) the space of all smooth complex- valued functions on M and the 
subspace of all functions with compact support, respectively. Assume that we are given a 
linear operator 

P : C^{M) C°^{M). 

If G is some chart in A4 (not necessarily connected) and k : G — > f/ its diffeomorphism 
onto an open set U C R", then let P be defined by the diagram 

C(f (G) — ^ G°°(G) 

K I Ik 

(note, in the upper row is the operator o P o iq^ where iq is the natural embedding 
iG : C^{G) ->■ C^{M) and ra is the natural restriction ra : C^{M) C°°{G); for 
brevity we denote this operator by the same letter P as the original operator). 

Definition 2.3. An operator P : C^{M) — > C°°{Ai) is called a pseudodifferential 
operator on Ai if for any chart diffeomorphism k : G U , the operator P defined above 
is a pseudodifferential operator on U. 

Lemma 2.4 (see, for example, Proposition 0.3.C of [?D]) // A and B are two 

pseudodifferential operators of order m andm' , respectively, then the composition G = AoB 
is a pseudodifferential operator of order m + m' with the symbol 

(2.5) c(x, ^ E ^^^(^' ^)D:b{x, 

^ — ' a! ^ 

a 

where a(a;,^) and b{x,^) are the symbols of A and B, respectively. In particular, the 
principal symbol of A o B is aQ{x,^)bo{x,^), where 09(2;,^) and bo{x,^) are the principal 
symbols of A andB, respectively. 
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Lemma 2.5. Let 



(2.6) 



A 



a-^"-- \ 



• a"-i^"-i 
V ••• a""/ 

&e a positive definite, real symmetric constant matrix. Let (t>{x') and h(x') be C°° Junctions 
of compact support in {n — l)-space. Then the problem 



u = 



on aM!^ 



(2.7) 

has a solution 

(2.8) u{x\xn)= I Ki{x' -y',Xn)Hy'W + [ K2{x' ~y',xn)h{y')dy\ 
where R" = {a; = (a;i, • • • ,a;„_i,x„) € R"|x„ > 0}, and 

n-l 



Ki{x',X,,) = (-1) 



(27ri)"-i 
+ (n — l)ixn 



|r,'|=l 



x' ■ rj' + ixr, 



„nn 



* X ■ rj + IXn 



a"" 1 



if2(s',2:n) = (-1) 



n-l (^-2)! 



(27ri)^ 



|r,'|=l 



X ■ 7] + IXr, 



dSri' ■ 



Here rj' — {rji, ■ ■ ■ ,77,1-1) and ds^r is the area element on the unit sphere \rj'\ = 1. 

Proof. Writing a; = (a;',a;„). Then the bi-Laplace operator P has characteristic form 
P{ri',T) — (X]j + a^^T^)^. It is easy to see that the roots of P{i]',t) with 

positive imaginary parts are T^{ri') = T2{rj') = i k=i "^g^n^ - Thus we have (see, 

Chapter I, §1 of [T) 



M+(r7',T) = 



n-l 



> a"" 



- 2i 



> a"" 



E 

j,k=l 



a^'^VjVk 



so that 

Ni{r,',T)=M+{T]',T)=T-2i 



\ j,k=l 



a^^riiVk 



n-l 

E 



\ 3,k=l 
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It follows from p. 635 of and the well-known residue theorem (see, for example, p. 150 
of [3]) that 



Ki{x',Xn) — { 1)" ^ (2,ri)" i|r)'| = l M+{rj',T)(x'-rj'+x„Ty 



—dr 



dSri' 



— ("-2)! r 

^ y (27I-1)" J|j7'|=i 



-dr 



dsjji 



(-1) 



n-l («-2)! 



1 iln'l^l 



1-n 



/ 1 \n-2 (»-3)! r 

(27ri)" J|ri'| = l 



/ lyi-l ("-2)! r 
V (27ri)"-i J|»)'l=l 



J 7 Af + (77' ,r ) (aj'-r/' +3:^1 r)"^-^ 



-dr 



X' •77' + «2;n-y/E",fcil 



ds„i , 



where 7 is a Jordan contour in Imr > enclosing all the points j=i "^^T^^' for all 



1 77' I = 1. Applying Theorem 2.1 of Jj, we obtain (|2^ . □ 

Let {Et} be the spectral resolution of pseudodifferential operator P, and let e{x,y,T) 
be the kernel of Et- This is an element of C°°{ft x fl) called the spectral function of P. 

The following Lemma will be used later. 

Lemma 2.6 (Homander's spectral function theorem, see, Theorem 5.1 of [T?], [H] 
or [32]) Let P be a non-negative pseudodifferential operator, acting on a C°° subdomain 
n of an n- dimensional C°° manifold. Let pq(x, S,) be the principal symbol of P, which is a 
real homogeneous polynomial of degree m on the cotangent bundle T*(D,). The measure dx 
defines a Lebesgue measure d^* in each fiber ofT*(fl); which is a vector space of dimension 
n. Then 



(2.9) 



r-"/'"e(x, X, t) - (2^)-" / = 0(t-i/") 



as T ^ 00, 



where = G T*in)\poix,0 < !}• 



3. The principal symbols 



3.1. Let {Ai,g) be a C°° Riemannian manifold, and let £7 be a bounded domain with 
C°° boundary in Ai . The "Neumann-to-Laplacian map" is the map 
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defined by the following problem: Let h e H^/^{dVl) and let u € H^{fl) be the solution of 



(3.1) 



A^-u = in n, 
u — on dil, 
P^^h on dn, 



we set Fh :— {—Agu)\g^. Multiplying p.ip by u, integrating the result over fl, and using 
Green's formula, we derive 

= / u{Alu)dx^ f \A,ufdx- f u^^4^ds+ f {A„u)^ds 
Jn ^ Jn Jan dv Jq^ du 

f f du 

= / \Agufdx+ / {Agu)—ds, 
Ja Jon <jv 



so that 

{Fh,h) = I {Fh)hds^ I \Agu\^dx>0, for any heH^^^{dn). 
Jan Jn 

This shows that F is a non- negative, self-adjoint, pseudodifferential operator on H^/^{dVl). 
We shall calculate the principal symbol of F. 

Lemma 3.1 Let A he a positive definition, real-valued constant matrix as in i2.6]) . Assume 
that 

Fo : C^(R"-i) ^ C°°{W~^) 

defined by the following problem: Let h £ {M.""^^) and let u G C°°(]R") be the solution 
of 



(3.2) 



u = 



we set Foh := - (E"fcii a^''-^ 
Po{x',r,')^2 



-,nn d 



on dW%. 

Then the principal symbol of Fq is 



n-l 



a^'^rjjTjk, 



V (x',/;') e M"-i X (M"-i \0). 



Proof. Writing x — {x' ,Xn), it follows from Lemma 2.5 that 
(3.3) u{x',Xn)^ I K2{x' - y' ,Xn)h{y')dy' 

where 



K,{x\x,,) = {-IT-'^^-^^ f 



x' ■ r( + ixrn 



dSrji ■ 
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(|3.3p shows that u{x) is uniquely determined by the data of h on its support set. p.2p can 
be rewritten as 



(3.4) 



Ik 

dxj dx^ 



u = 

■J Qnn du 



on dRl, 
on 9R" . 



Taking the Fourier transform to p.4p with respect to xi, • • • , Xn-i, we have 



(3.5) 



We denote 



n-l \^ 



u = 



/i(v) 



on dR^, 
on aM'l. 



Then, the general solution of p.5p has the form: 



(3.6) uiT]', Xn) = Ciel^'l"" + C2e-I«'l"" + C3X„el«'l"" + C4a;„e-I«'l"", 

where Ci, C2, C3, C4 are arbitrary functions in r]' . From the boundary conditions of p.Sp 
it follows that 

m(77',x„) = C2(-el«'l^"+e-l«'l^"+2a;„|C'|e-l«'l"") 

+C3(x„el«'l-" - a;„e-l«'l-") + ^x„e-l«'l-. 



Therefore 

u(a;) 



(27r)^ 



from which we have 

d^u ] 



dxjdxk (27r)"-i jRn-i 



e'^"'''''>(-?7,?7fe) fc2(-el«'l"" +e-l«'l^" +2x„|e'|e-l«'l"") 



(27r 



\n — l 



+2x„|C?e-l«'l^") + C3 (2|e'|el«>" + x„|r pel«>" + 2|r|e-l«>" - pe"!?' 
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Then 



dxjdxk 



dxl 



a: ,1=0 



dxl 



{2n) 



x„=Q 

n Hi) 



In order to take a bounded solution of the equation p.Sp . we may let C2 ~ C-z — Q. Hence 



.Ik, 



dxjdxk 



'1^ 



n HV') 



1 



(27r)»-i V-i 
1 



-2 



I.e. 



(3.7) 



{2n) 



Foh 



-2, 



E ^^''VjVk I e 
\ j,fc=i 



(27r) 



This shows that the principal symbol of the pseudodifferential operator Fq is 2 y fc=i ^.-'''VjVk- 
□ 

Theorem 3.2 Let (A^,.g) 6e an n- dimensional Riemannian manifold, and let he 
a hounded domain with C°° houndary. Assume that the pseudodifferential operator F is 
defined as before. Then for any coordinate chart k : dfl D U ^ U*^ C R"^^ there is a 
pseudodifferential operator A G '^^{U'^) such that for every h G H^/'^{dQ.) we have 

K,{Fh) ^ A{K,h) e C'^iW), 

where is the linear tangent mapping of k. A is elliptic and has a coordinate invariant 
positively 1-homogeneous principal symhol pq G C°°{T*U'^ \ 0) = C°°{U'^ x (R""^ \ 0)) 
given hy 



J2 a'^mm, V(a;',r/)GT*{/^\0. 



Proof. It is well-known (see, for example, [24|) that there is a T > and a neighborhood 
G C of the boundary dfl together with a diffcomorphism i/) : G — > dfl x [0, T) such that 

i) tplq) — {q, 0) for every q G dfl, 

ii) The unique geodesic normal to dQ (with the unit-speed ^/g™' with respect to g) 
starting in g G i9il is given by 

; } • 

Moreover, ip is unique with i) and ii) and has the following additional properties: Let 
K : on D U -i' C R""i be any coordinate chart on dn and k : M ^ IJ x [0,T) 

be its extension via ip. 
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1) The metric g has on x [0,T) the form 

n-l 

(3.8) K*.g = ^ [gjk dxj (g) dxk) + ff„„dx„ (g) (ia;„, 

i,fc=i 

2) The Laplace-Behrami operator Ag can in x [0,T) be written as 



5 



1 5(Vl5l5"") 



k is cahed a boundary normal coordinate chart and its coordinates xi, ■ ■ ■ , x„_i, a;„ bound- 
ary normal coordinates. G is said to be a tubular neighborhood of d^l. Therefore, for given 
e > and every q G dQ, we let {Gq^^,^) is a boundary normal coordinates chart (Note 
that k{Gq^e) = Uq ,: X [0,T)) such that diam(G5,e) < e. Then there is a partition of unity 
subordinate to the open cover {Gq^e H dn\q € Oil}, i.e., a collection of real- valued C°° 
functions 4'i on dfl satisfying the following conditions: 

(a) The supports of the (f)i are compact and locally finite; 

(b) The support of is completely contained in Ga for some a; 

(c) The (j>i sum to one at each point of dfl: 



1. 



Since h — h(f>i, we may assume that the support set of h is contained in some small 
neighborhood Gq^^ H dil for some q G dil. (Let us point out that (^(G'g^e)) n {x„ = 0} = 
Uq ,: C 9R" ). It is clear that we can always choose a fine cover {Gg^e n di}} of dfl, so that 
in addition to p.Sp we have 



(3.9) 



\{g'\x\x„)-g^\0))T^,T^k\<eJ2 



{{x',x„) e k{Gq^,)) 



for any given e > 0, all real 771, •• • and all k{Gq^^). The finer cover does not influence 
the fact, which is obviously true in the original cover, that p.9p holds. 

Let u be the solution of p.ip . For any fixed eo > 0, in the local coordinates, the u 
satisfies 



(3.10) 



E 



J,fc=l 



d 



+ 



+ 9 
a(Vlg|9"") 

dXn 



nn d 
dx'i, 



d 



in K(G,,eJ, 

on k{Gq^eo) n I 
on k^Gq^eo) ' 



By the regularity of elliptic equations, we get u G G*'"(K(Gq_e)). We define the operator 



<3,eo ■ C'5"((K(Gg,co)) n 



^ C-((^(Gg,eJ) n 9R!^) by Fq,,,h (-Au)!^^^^^ ^^^^^ 



for any h G G^((K(Gq,eo)) n 9M"). It is easy to check that Fq^eo is a pseudodifferential 

operator on the open set (^(Gq^eo)) dW^- We denote its principal symbol as p^Q°\x', rj'). 
Next, noticing that (see, for example, Theorem 1.4.4 of |T5]) 

dg^'' 



dxi 



(0) = for aU l<j,kj< n, 
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we get 



92 



dxjdxk 



g""(0) 



where is the boundary normal coordinates of g G dfl. Let v be the solution of the problem 

2 

w = 



(3.11) 
Since 



(E"^ii5^'=(0)g^+5""(0)^)% = in R'^. 



V9^m = h 



on 
on 



we get that v — u satisfies 



"9(0) 

V — u = 



Al)u in k{Gg,,), 

in K(Gg^e), 
on (K(Gq,,)) n ^R'J: 



V9^)^^^iV9^)~V9^))^ on (^(G,,,))naM;. 

It follows from i^-estimates of elliptic equations (see, for example, Theorem 15.3 of [1]) 
that 



\/g""(a:)--v/g""(o) a-u 



+ Ik - w||lp(S(G,,.)) , 



where the constant Ci is independent of w — u, W'"P{K{Gq^e)) is the Sobolev space with 
p > n, and < e < eg. From 



lim \\v~u\\lp^^g,.,)) =0, Imi |1 (A^ - A2(oj)u||ip(s(G,,,)) =0 



and 



lim 

e-)-0 



\/<?""(a:) - v/5""(0) 9m 



'(0) 



0, 



we obtain 



Combining this and the Sobolev imbedding theorem, we find 

|jw-u||c3.-(;j(G,,,/2)) =0- 

Furthermore, applying the Schauder estimates (see, for example, the proof of Theorem 7.2 
in [T]), we have that for any e < eo, 

\\V - w||c^-(S(G,,,/4)) < ^2 (||(A2 - A2(o))7,||c.(s(G,.,/.)) 



c3'°(S(G,,,/2))naK!; 



+ I1^'-w||c-(k(g,,„)) , 



where the constant C2 is independent of 1; — u. By 
1™II(A'-AL))« 



g(0);"IIC°(S(G,,,/2)) 
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and 



lim 

e->0 



^.g«"(a;) - 1/5"" (0) du 



5""(0) aa;„ 



0, 



C3'°(K(G,,,)naR:j:) 



have 



which implies 



l™ll^-^llci-(«(G,„/4)) =0' 



lim II Agu - Ag(o)w||c2,»(s(G,„/4)) = 0- 
Applying the trace theorem, we obtain 



lim 



I ((^s'")|c(K(G,,,/4)naR!^) (^9(0)")lG(K(G,../4)naR!^)) 



= 0, 



I.e., 
(3.12) 



{x'=0 (x',o)eR"-i 



It follows from Lemma 3.1 that has the principal symbol 2^ 9-'''iO)VjVk, where 

Fq : C^ikiG,^,„)ndRl) ^ C°°(^(G,,eo)n9M'^) defined by Foh = (A«)| _(^^^^^^)^gjj„ , and 
V is the solution of p. lip . From p.l2p . we get that the principal symbol of Fg,eo is also 
p'o°\^i^') — ^A/X^^feii 5"''' (0)??j?7fc- Since q is an arbitrary point at dfl, it follows that the 



principal symbol of F is 2^J2jk=i9''''(^')VjVk in the local boundary normal coordinate 
system (x') of dH.. □ 

3.2. We define another pseudodifferential operator Q : H^/^{dVl) H~^/^{dfl) as 
follows: Let (j) e H^^^{dn), and let u e H'^{9) be the solution of 



(3.13) 



we set 8/1 



U = (f) 
fi = 



in Q, 
on 90, 
on dfl. 



By Green's formula, we have 



ds= \Agu\^dx > 0, for any h e H^^^{dn), 
Jn 



{eh,h) = / {eh)hds ^ I 

Jdn Jan \ 

which implies that 8 is a non- negative, self-adjoint pseudodifferential operator from H^/^{dVL) 
to H-^/^dfl). 

Lemma 3.3 Let A be a positive definition, real-valued constant matrix as in i2.6\) . Assume 
that 

80 : C^iW-'^) C°°(M""i) 

defined by the following problem: Let (j> £ (W^^^) and let u £ C°°(R") be the solution 
of 

2 



(3.14) 



^nn d 
dxt 



' U = (f) 

OXn 



u = Q 



on 
on 
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we set Qo(j) := V a"" ^^^""^ \dR" ' -^^^^ ^'^^ principal symbol of is 



n-l 



V (a;',ry') G M"-^ x (M"-i \ 0). 



Proof. Similar to Lemma 3.2, it follows from Lemma 2.5 that 
(3.15) u{x',Xn)^ J Ki{x' - y' ,Xn)<j){y')dy' , 

where 

^n-l (^-2)! 



Ki{x',Xn) = (-1) 



+ {n — l)iX' 



x' ■ rj' + ixn. 



> a"" 



Taking the Fourier transform for p.l4p with respect to xi, • • • , a;„_i, we have 

2a'' 



dSr, 



(3.16) 



u = (j){r]') 
#^ =0 



+ (EI fc=i a^Sj'yfc) ^ = in R!^, 

on d. 



on dWl. 



Then, the general solution of p. 161) has the form: 

(3.17) u(i, Xn) = Ciel«'l"" + C2e-I«'l"" + C^x.J^'^''- + C4a;„e-I«'l"", 



where |^'| :— \/ ^^''° , and Ci, C2, C3, C4 are arbitrary functions in 77'. In order to 

obtain a bounded solution of (I3.16p . we put Ci ^ Cz ~ 0, so that we find by the boundary 
conditions of (jXTCI that C2 = 0(?/), C4 = 0(r?')IC'l- That is, 

zl(7;',x„) = <^(r;')e-l«'l^"(l + x„|^'|). 

e*(^',';'>0(^')e-l«'l=^" (1 + x„|C'|) rf??'. 



Thus 



u{x) 



i2n) 



-1 



Since 



1 



dxjdxkdxn 



dxl 



(27r)"-i 
1 



e*<-'"'>0(r;')e-l«'l-" (2|C? - x„|Cr) d^?', 
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it follows that 

eo<? 



ox, 1 



dxjdxk 



dxl 



1 



1 



*<^'^'''>0(ry')d?7'. 



This shows that the principal symbol of the pseudodifferential operator Qo on dMJf is 



VjVk 



□ 



Theorem 3.4 Lef {M,g) be an n- dimensional C°° Riemannian manifold, and let Q be 
a bounded domain with C°° boundary. Assume that the pseudodifferential operator is 
defined as before. Then for any coordinate chart k : dfl D C/ — >■ U'^ C K""'^ there is a 
pseudodifferential operator T e '^^{U'^) such that for every (j) e H^^^{dVL) we have 

T is elliptic and has a coordinate invariant positively S-homogeneous principal symbol po G 
C°°{T*U^ \ 0) = C°°([/'^ X (M"-i \ 0)) given by 

3/2 



Po 



ix',rj') = 2 ^ g^'^Wk 



V ix'.r]') e r*c/'^\o. 



Proof. The proof is similar to Theorem 3.2. 



4. Proofs of main results 



Proof of Theorem 1.1. Let F : H^/^{dn) H-^/'^{dn) is defined as in §3.1. It follows 
from the discussion in §3.1 that is a self-adjoint, elliptic, pseudodifferential operator on 

whose principal symbol is 2^^"^^^ g^^{x')rijrik, where x' is the local boundary 
normal coordinate on dfl. We define the operator by 



Zefix') 



1 



fix') for all / e H^^^idn) and x' e dfl, 



^g{x') + e, 

where e > is a sufficiently small constant. Applying Lemma 2.4, we obtain that the 
operator Q, = Z,oF : H^/^idfl) ^ H-^/^id^) defined by Q,h = ((^(^^ (-Agu)) 

is a pseudodifferential operator with the principal symbol — - — g(x')+e ' ^h^'"^ 

is the solution of p.ip . It is easily seen that the operator has the same eigenvalues 
Xk (e) and corresponding normalized eigenfunctions the following biharmonic 
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Steklov eigenvalue problem: 

Aluk = in n, 

Uk — on dfl, 

AgUk + {Xk{e)){g{x') + e)^ = on dn. 

Let Er be the spectral resolution of Q^, and let e{x',y',T) be the kernel of Et (here 

eix',y',T) = Ea.<x + ^)iQiy') + e) ^^)- K follows from Lemma 2.6 

(see, Theorem 1.1 of [II], or [12]) that 

(4.1) r-("-i)e(a;',x',T) - (27r)-("-i) / d^* = 0(t-1) as r ^ +c5o, 

where B^' = {rj' G T*,{dn)\po{x' ,1]') < 1}, po{x',ri') denotes the principal symbol of Q^- 
By A^{t) — e(x', x' , T)dx' and Pq{x' , -q') ~ ^^^'g^a^+e ''''^ ^^ we obtain that 

Mr) = j,r^Al dx' [ dc]r-~' 

+0(t""^) as t ^ +00. 

For any fixed local boundary normal coordinate x' e k(dQ), since (ti — 1) x (n — 1) matrix 
g' = {g^^{x')) is positive definite, there exists an (n — 1) x (n— 1) matrix C(a;') = {cjk{x')) 
such that *C{x')g' {x')C{x') = (Sjk), where Sj/. is the Kronecker delta. Note that d^* = 
■\/|.g'(x')| dCi • • • dCn-ij in each fiber of T*(9f2), which is a vector space of dimension n~l. 
With the change of variables rjj — X]fc=i Cjk{x')Ck, we obtain 

/ dC 



, , |detC(x')| v^7MrfCi---rfCn-i) 

{(Ci,-,C,.-i)eR"-V?i'+-+^"-i<^^^> 

, d(i-- ■ dCn~l 

{(Ci,-,C.-i)eK"-i|^Cf+-+CS_i<^^i+i} 

;^.(x') + e^"~' 



here we have used the fact that |detC(a;')|-\/|g(a;')| = 1, and where a;„_i is the volume of 
the unit ball of R"-^ Therefore 

1 ,. ri-l f ( + J I , r^^ n-2\ 



Letting e — > 0, we obtain 

^(r) = ^^^^c.„_ir"-i / 'dx' + 0(r"-2) ^s r ^ +oo, 



an 



that is, 



Mr) = ^;;^^„_i / igis)r-'ds + 0{t-') as r ^ +«) 



□ 
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Proof of Theorem 1.2. Let Re : H^/^{dn) H-'^/'^{dVt) be defined as follows: For 

any </> g H^/'^{d^), we put Re<p — ((^q^ '^^aj"^ ) I an' where v satisfies 

Alv = in n, 
u — (j) on (9i7, 
= on dn, 

and e > is a sufficiently small constant. Clearly, is a self-adjoint, elliptic, non-negative 
pseudodifferential operator of order 3. By Lemmas 2.4 and Theorem 3.4, we get that 

the principal symbol of R^ has the form (^ipjp ( fc=i 5"''^'7j^fe ) ' where ?/ e R"^^. 
It is easily seen that the operator R^ has the same eigenvalues /ifc(c) and corresponding 
normalized eigenfunctions Vk on Sfi as the following biharmonic Steklov eigenvalue problem: 

A^wfc = in f], 

1^ = on 
%^-/i^(6)(,-|-e)3^;fe = on 9a 
It follows from Lemma 2.6 (Theorem 1.1 of [13], also see [12] or [5]) that 

(4.2) r-("-i)/3g(^/^^/^^) _ (27r)-("-i) / d^* = 0(t1/3) as r ^ +c5o. 



where Bj;' = e T^,(9ri)|po(2:','7') < I}, Pq{x',v') denotes the principal symbol of F^. 
Since Be[T) = J^j^ e[x ,x , r-'jdx and po(a; , ?? j = -^^ (g(x')+e)3 ' '^^^'^ 

B,{t) = I _^ ( / da;' / 3„ ) t"~' 

(27r)« J2(e(x')+e)-3(y'5:;-i^5J'--(:E')r,,„fe) <i / 

-|-0(t""^) as t ^ -fcx), 

where d^* = \/W{^') \ dQi ■ ■ ■dC,n-i- With the change of variables r^j = X]fc=i Cjk{x')(k, 
where (n — 1) x (n — 1) matrix C{x') = {cjk{x')) satisfies *C(x')(?'(x')C(x') = {6jk), we 
obtain 



dCl • • ■ rfCn-l 



which implies 



n-l 



n-1 



Letting e — > 0, we obtain 



I.e., 



A{r) = / e"-i(s)ds + 0(r"-2) as r ^ +c^. 
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□ 



5. The asymptotic formulas are sharp 



Hormander (see [14] or [32]) proved that, for a pseudodifferential operator of order in 
with principal symbol po{x, f ), 



(5.1) 



e{x,x,T) - (27r)-" / dC 

Po(a;,?)<T 



<cii + \T\y- 



uniformly in compact subsets of Ai, where C is independent of x,t. Applying this result 
to our cases, we immediately obtain 

^ri-l 



(5.2) 



(5.3) 



Air) 



(4^)"-' Jdn 



<C(1 + |t| 



\n-2 



t^ri-lT" 



(^^)"-l JoD 



<C(l + |r| 



In this section, we shall show that ()5.2p and (|5.3p cannot further be improved. More 
precisely, we shall show by two counterexamples (letting respectively be the unit ball of 
M" and the unit disk of of IR^,and g= I) that the asymptotic formulas (11.81) and (II. 9p are 
the "best possible". 

First we give some well-known facts concerning spherical harmonics (See e.g. Miiller 
[5S]). When il — B we may explicitly determine all the biharmonic Steklov eigenvalues 
of p.ip . In fact, for each integer m > 0, let 7'm(R") denote the set of homogeneous 
polynomials of degree m in n variables, i.e., the set of functions u of the form 

u{x) = ^ aaX°' for x G M", 

\cx\—m 

with coefficients aa G C. A solid spherical harmonic of degree m is an element of the 
subspace 



(5.4) 

Let 

(5.5) 



■H™(R") = {m e ■P™(]R")|Au = on K"}. 



N{n,m) = dim-H,„(]R") for n > 1 and m > 0. 



Note that = Hq is just the space of constant functions, and Vi = Hi is just the space 
of homogeneous linear functions, so 

(5.6) N{n,0)^l andiV(n, l) = n for n > 1. 

It follows from p. 251-252 of 22 that 

1 if TO = or 1, 
if TO > 2, 



iV(l,TO) = 

7V(2,to) = 



1 if TO = 0, 

2 if TO > 1, 



SHARP WEYL-TYPE FORMULAS 



19 



and 



N{n, m) 



2m + n — 2 f m + n — 3 



n-2 



71 — 3 



for n > 3 and to > 0. 



The following Lemma was obtained by Ferrero, Gazzola and Weth (see, Theorem 1.3 of 

my- 

Lemma 5.1. If n > 2 and Vl — B, then for all m — 0, 1, 2, 3, • • • ; 

(i) the eigenvalues of are Am — n + 2m; 
(a) the multiplicity of Xm equals N{n,m); 

(Hi) for all tpm € Hmi^""), the function (j)m{x) '■= (1 — |xp)'0m(a^) is an eigenfunction 
corresponding to A,„. 

Now, let < Ai < A2 < • • • < Afc < • • • be all biharmonic Steklov eigenvalues for the 
ball B. From the above lemma and the formula 

^ ( a+ j -I \ _ (a + l)(a + 2)---(a + TO) _ 
^ I J ; m! 



we get that for n > 2, 



fc=0 



^ (7i + 2fc-2)(n + fc-3)! _ ^ 

k=l ^ ' k=l 

n + k-2 \ (n-2)! " 



(n + k -2)1 (n + fc-3)! 



= 1 + E 

m 



/ (n-2)! 

k=l ^ / \ J 1.^2 



n + fc — 3 
fc- 1 



(n-2)!fc! (n- 2)!(/c- 1)! 



k=l 



(n - 1) + - 1 
k 



1 + E 



k=2 



(n- 1) + (A: - 1) - 1 
k-1 



2 + E 



(n- l) + k-l 
k 



-E 

k=i ^ ' j=i 

n(n + 1) • • • (n — 1 + to) 



(n-l)+j-l 
j 



(n - 1) + to - 1 



1 



n(n + 1) • • • (n — 1 + to) 



n + TO — 2 

TO 



= 2 
= 2 



n(n + 1) • • • (n — 1 + to) 



TO 

n + TO — 1 

TO 



n + TO — 2 

TO 



rt + TO — 2\ 2^ n + TO — 1 
TO / V 71 — 1 



n + TO — 2 

77-2 



By applying the formula 



p+l 
r 



P 
r - 1 
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and m — ^ — ^, we obtain 



n + m — 1 
n- 1 



71 + m — 2 
n- 1 



{2m + n - l)(m + n - 2)(m + n - 3) • • • (m + 1) 

{n-iy. 

1 _ _ 

(A„ - 1)(A„ + n- 4)(A™ + n - 6) • • • (A„ - n + 2). 



2"-2 (n - 1)! 
In view of ^^^^^ = r(z)r(z + i), we get 
1 2"7r"-i 



2"-2(n-l)! (47r)"-i(n-l)! (4^)„-i £ir(n+i) 



1 



(47r)"-ir(f + i)r(f + 1) (4^)"-i r(2i2^ + i) r(f + i) 

(47r)"-i (47r)"-i^ ^ " 



Hence 
A(A™) 



(47r) 



(47r) 



^(vol(9S))(A„ - 1)(A„ + n - 4)(A„ + n - 6) • • • (A,„ - n + 2) 
^(vol(9i3)) [A(;r 1) + (1 - n)A:r2 + 4)(„ _ g) • • • (-n + 2) 



Since 1 — n 7^ 0, this shows that the formula (jl.Sp is sharp. 

Similar to the above yl(Am), we can also give a counter-example to show that the re- 
mainder term estimate in the asymptotic formula (I1.9P is sharp. Let _B C R'^ be the unit 
disk. If m > 1 , then the functions 



(5.7) 



"ip m,i{i', 9) — r"^ cos mO and iprn.2{'r', 9) — r"^ sin m9 



form an orthogonal basis of harmonic function in the space "HmlR^)- Let us consider the 
following Neumann boundary value problem: 



(5.8) 



1pm, j in B, 
= on dB, 



where m = 1, 2, 3, • • • ; j = 1, 2. We claim that the above solutions satisfy 

in B, 



(5.9) 
where 



^""-i = on dB, 

1 a(Au„.,j) 



9i/ 



on dB. 



JOB [ df 



ds 



JdB \ du J 



ds 



2m2(m + 1) 



In fact, for the Dirichlet problem 



Aum,] = ipm,j{r,9) in B, 

^m,j = 2m(m+l) V'm,j(l,6>) On DB, 



m ^ 1,2,3,- ■■; j = 1,2, 
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from the formula of the solution to the Dirichlet boundary value problem for the Poisson 
equation in the unit disk, we have 



1 



Urn,jip,a) ^ - — / rdr In 



47r 
1 



27r 



1 + p^r^ — 2rpcos(Q! - 



+ — 2rpcos(a 



2tt Jq 1 + p2 _ 2/9Cos(a - 6) \2m{m + 1) 



dO. 



Then 



dumjip, a) I 



1 p2tt 



r(l - r2) 



dp 'P-^ 2tt Jq Jq 1 + — 2r cos(a — 



-ilJm,j{r, 9)drd9 



1 

2^ 



d 



2-iT 



dp 1 + p2 — 2pcos(Q; — 6*) \lm(rn + 1) 
It is well-known that 

l-r2 



d6l 



7. -, rr = 1 + 2 r' cos - 

1 + r2 - 2rcos/(a - 0) ^ ^ 



Therefore 
dp 



lp=i 



1 ^1 ^27r / oo \ 

-y y f r + 2^r'+icosZ(a-6')j Vmj(r,6')drd6l 

1 r /-ZTr / oo \ 



1 

'2^ 
1 

'2^ 



1 ^Q-27r 



J a 



/a-2ir / oo 
2^;p'-icosZi 



2'm{m + 1) 
2 ^ r'+^ cos Itj Tpm.j {r, a - t)dr dt 



p=l 



1=1 



1=1 

dujj 



2m(m + 1) 



■dt 



By (|5.7p and a simple calculation, we get g™'^ |p-i = 0- Combining this and the result 



where 



in B, 

-ar~ - Vm,jip7n,j on dB 
1 /B|V(V'™j)l''^a: 1 



m=l,2,3,-- - ; j = l,2. 



we show the desired claim. It follows that 



— , TO = 1,2,3, j = 1,2, 



/i^, = 2to2(to + 1), TO = 1, 2, 3, • • • ; j = 1, 2, 



Thus, for 71 = 2, 



m ^ 
B{Pm.j) = <A'mj} = $]^(2,fc) = 2(TO+-) 

A:=0 

^ c.i(voi(ai?)) / 1 

27r V 2 
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By 



as m — > +00, 



we have 




as TO — )• +00, j = 1, 2. 



Since 5 7^ 0, this shows that asymptotic formula 



(|1.9I) cannot be improved on unit disk of 
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